Quantum destruction of stiffness in diluted antiferromagnets and superconductors 
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The reduction of 2D superconducting or antiferromagnetic order by random dilution is studied as a 
model for the 2D diluted Heisenberg antiferromagnet (DHAF) La2Cui_ p (Zn,Mg) p 04 and randomly 
inhomogeneous 2D suerconductors. We show in simplified models that long-range order can persist 
at the percolation threshold despite the presence of disordered one-dimensional segments, contrary 
to the classical case. When long-range order persists to the percolation threshold, charging effects 
(in the superconductor) or frustrating interactions (in the antiferromagnet) can dramatically modify 
the stiffness of the order. This quantum destruction of stiffness is used to model neutron scattering 
data on La2Cui_ p (Zn,Mg) p 04. In a certain simplified model, there is a sharp stiffness transition 
between "stiff" and "floppy" ordered phases. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Cx, 75.40.Mg 



Randomly diluted superconductors and antiferromag- 
nets involve a combination of classical percolation physics 
with the quantum physics underlying superconductivity 
and antiferromagnetism. Percolation is perhaps the sim- 
plest transition that can occur in a disordered system: 
random dilution of sites (or bonds) on a lattice induces 
a transition between a phase with one infinite nearest- 
neighbor connected cluster of occupied sites (bonds), and 
a phase with only disconnected finite clusters. Recent ex- 
periments on diluted 2D antiferromagnets and inhomoge- 
neous superconductors require a theory of how quantum- 
mechanical effects modify the percolation transition in 
these systems. This question is also of practical impor- 
tance for field-effect devices in which a thin film is tuned 
through the superconducting transition []J. 

At the percolation threshold, p c , the infinite connected 
cluster is on the verge of being cut into disconnected fi- 
nite clusters (Fig|TJ. For the connected cluster to have 
long-range order (LRO), the quasi-2D "blobs" must cor- 
relate across quasi-lD "links." This occurs easily in some 
other diluted magnets 0, and inhomogeneous super- 
conductors Q , since their degrees of freedom have LRO 
even in a ID chain. The two cases considered here, the 
s = 1/2 Heisenberg antiferromagnet and the O(N) quan- 
tum rotor, both have quantum degrees of freedom that 
order in 2D but not in ID. Then the question of whether 
LRO survives to p Cl when the infinite cluster is a fractal 
object of dimension || w 1.896 is unanswered. 

The effect of quantum fluctuations is strongest on the 
quasi-lD chains. For classical models on the percolation 
cluster at T > 0, the fact that these chains are disor- 
dered (have a finite correlation length £) implies that 
the cluster has no LRO. Similar arguments have been 
made for the quantum cases discussed here Q . The first 
part of this paper shows that the existence of arbitrar- 
ily long ID segments at p c , and the fact that a spin or 
rotor in the middle of such a segment has strong quan- 
tum fluctuations, does not prevent LRO for the cluster: 




FIG. 1: The backbone of the incipient infinite cluster (the 
backbone is the portion that carries current from one end of 
the sample to another) showing 2D "blobs" and ID "links" . 

the 2D blobs can order through disordered ID links, in 
a manner that is impossible for classical models. These 
one-dimensional segments cause difficulty for spin-wave 
calculations because, as seen below, the spins in the mid- 
dle of such a segment fluctuate strongly. An additional 
motivation for studying the effect of these ID segments 
is that their effects are essentially unobservable in cur- 
rent QMC studies, as a typical realization of e.g. 10 5 
spins will contain no ID segments of length longer than 
8. A different way of connecting ID physics to diluted 
antiferromagnets is discussed in |6J. 

Randomly diluted quantum degrees of freedom ap- 
pear in two well-known nanoscale inhomogeneous materi- 
als. La2Cui_p(Zn,Mg) p 04 is obtained by adding "static 
holes" (i.e., removing spins) at random in a quasi-2D 
antiferromagnet. For small hole densities pCp Cl neu- 
tron scattering measurements |j| agree well with quan- 
tum Monte Carlo 0, H, and spin- wave calculations 
using the s = | site-diluted Heisenberg antiferromagnet 
(DHAF) model 

H = J^piPjSi-Sj, (1) 

where the sites form a square lattice, and Pi,Pj equal 
with probability p and 1 with probability 1 — p. (Vajk 
et.al. argue that other interactions besides those 
in (JJJ may not be neglected close to threshold. We con- 



sider some such interactions below.) 

The second type of material, inhomogeneous high- 
temperature superconductors like E^S^CaC^Og+s 
(BSCCO), comes from doping mobile holes into a quasi- 
2D antiferromagnet. Scanning tunneling spectroscopy of 
BSCCO surfaces show grains of size ~ 3 nm that are ei- 
ther superconducting or insulating |l2j . The grains are 
larger than the coherence length w 1 nm, so it is rea- 
sonable to assume that they contain Cooper pairs, and 
that collectively they resemble a bond-diluted Josephson- 
junction (JJ) model 



Hjj = J2 E c(ni 
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Here the bond variables pij have the same distribution 
as pi in JU . The above is not expected to be as accurate 
a model of the microscopies as the diluted Heisenberg 
model Q is for La2Cui_ p (Zn,Mg) p 04, since the micro- 
scopic origin of the disorder in BSCCO is unknown. 

The model (01 is in the universality class of the 0(2) ro- 
tor model, and for Ec ^ 0, charging effects in the grains 
prevent LRO in a ID chain 0] . The local charge rii does 
not commute with the superconducting phase (f>f. charg- 
ing tends to fix the local number, which is conjugate to 
the local phase. In the zero-charging-energy limit Q, 
Hjj orders in ID at T = 0. In the Heisenberg spin-half 
case, and for small nonzero Ec/Ej in the Josephson- 
junction case, the ID model is critical (has power-law 
correlations); for large Ec/Ej or Heisenberg spin-one, 
the ID model is short-ranged. 

One of our main conclusions for both the Heisenberg 
antiferromagnet and O(N) rotor is that order at p c is 
allowed despite the existence of ID links: in a simpli- 
fied limit we find a "renormalized classical" phase that 
provides an adequate description of some previous nu- 
merical results. We find that possible ordered states at 
p c must have extremely low stiffness relative to the undi- 
luted 2D case. At T = 0, order and stiffness are inde- 
pendent quantities, but mix as T increases: low stiffness 
leads to a reduced T > correlation length A. Such or- 
der with low stiffness can be proven to occur if there are 
classical magnetic moments or superconducting phases 
(see ©) coupled by quantum ID links. Superconducting 
phases with small superfluid density have previously been 
proposed, e.g., the "gossamer superconductor" of [l4| . 
Our picture differs in that the low superfluid density re- 
sults from randomness on scales larger than the coherence 
length [l2l llf^ , rather than from a uniform theory. 

We first calculate correlations in a toy model JSJ to 
show explicitly how quantum disordered or critical T = 
ID systems like the spin-half chain are fundamentally dif- 
ferent from classical disordered T > ID systems. Even 
though both may be disordered in ID, the quantum sys- 
tems can order through disordered ID regions (FigGJ 
while similar classical systems cannot. This simple state- 
ment underlies the possible existence of order at p c :. We 
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FIG. 2: Two examples of how quantum systems can order 
through a disordered ID segment. The first case is an AF 
spin-half chain with two weakly connected spins at the end; 
the second case is two bulk superconductors (Ec = 0) con- 
nected by a chain of grains with nonzero Ec- 

next solve a rather simplified limit of J2J) that shows a 
stiffness transition between two ordered phases: a "stiff" 
renormalized classical phase, and a "floppy" phase dom- 
inated by fluctuations on the ID links. In a less sim- 
plified but still approximate model, the "stiff phase" be- 
comes just a renormalized classical phase like in the 2D 
Heisenberg antiferromagnet (AF) except that the 

renormalized problem is not the 2D classical AF but the 
classical AF at p c . Finally, we compare experiment and 
simulations of to this renormalized classical theory. 

First consider a chain of spins with uniform coupling 
J, plus one spin at each end attached by a coupling J': 



JV-l 

Hi = 



Bi+l + J'(S ■ Si + Sat • Sjv+l). 



(3) 



Here, Sj are s = | Heisenberg spins coupled antiferro- 
magnetically, J ^> J' > 0. For very small J', the state of 
the internal N spins is nearly undisturbed, and in par- 
ticular si and s^v are only weakly correlated with each 
other However, the two spins at the ends So and sn + i 
can be made to form a perfect singlet with each other, 
for N even: limj7/j_ > o+ s • s^v+i = —3/4. 

To understand this result, note that at J' = 0, the 
four states of the end spins, which can be classified into 
a singlet and triplet, are degenerate. Perturbation the- 
ory in J' preserves rotational symmetry and, for J' <C J, 
just splits the singlet and triplet. The two spins at the 
ends correlate perfectly with each other, even though the 
internal spins are weakly correlated. Diagonalization, us- 
ing thick-restarted Lanczos [13, up to N = 20 confirms 
the second-order perturbation theory argument. At small 
enough J'/ J there is an arbitrarily strong singlet corre- 
lation 
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but with a vanishing energy splitting between singlet and 
triplet: 

( J'\ 2 

^triplet — ^singlet ~ /a(iV) J I — J . (5) 
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Here /i and are dimensionless functions of the chain 
length that incorporate the matrix elements and energy 
denominators of the unperturbed chain. 

At p ci the Josephson junction array J2J), in a certain 
approximation, also exhibits long-range order. Consider 
(for now) neglecting all fluctuations within 2D blobs. The 
scaling properties of this model can be found exactly; 
later, fluctuations within the blobs will be partially re- 
stored. In this approximation, the phases of the blobs, 
{<pi}, commute with H, so choose a basis of simultaneous 
eigenstates of (01 and <fii. In this basis, the energy 
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depends on exchange constants Ji > chosen from a dis- 
tribution P{J) which we calculate below using properties 
of the ID links (see 0). Regardless of the distribution, 
we have long-range order, since Ji > and <pi are classi- 
cal. We calculate the stiffness of this order below using 
the J -y behavior of P(J). 

Consider a single ID link of L sites bounded by two 
superconducting blobs of phase <fii and <j>2- After the 
quantum-classical mapping, we obtain a 2D XY model 
with fixed boundary conditions |l3l | . The lowest energy is 
obtained if <pi — 4>2, but the energy cost to create a phase 
difference depends on the coupling K = Eg J E.j in the 2D 
XY model. This problem is discussed in |l8| : the stiffness 
per site goes to a finite limit in the algebraic (ALRO) 
phase of the 2D XY model but falls off exponentially 
with L in the short-ranged phase (SRO). For small phase 
difference between boundaries Atfi = 4>i — 4>2, 
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The exchange strength is J(L) = 2E{A(j>) / '(A</>) 2 . 

The distribution of exchange strengths between blobs, 
P(J), is chosen to reflect the geometry of the percolation 
cluster 0, 0|. The total number of ID links between A 
and B goes as N ~ i? 3 / 4 and the fraction of links P(L) 
of length L falls off exponentially with L: P(L) ~ e~ aL . 

The probability distribution P(J) of exchange 
strengths is a sum of 6- functions: 



P(J) 
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P(L)S(J - J(L)) 



(8) 



The peaks become closely spaced as J — > 0. Using JJJ) to 
change variables, these properties imply that as J — > 0, 
the fraction of links with exchange J is 



p( T n s / fci/J 2 e- fclQ / J in ALRO phase 
1 j \ (J/ki)^- 1 in SRO phase 



(9) 



Any given link must be shorter than the total number of 
links N ~ i? 3 ' 4 , so the distribution © is cut off at 



T _ ( ki/N in ALRO phase 
1 k 2 e~ N /< in SRO phase 



(10) 



Now consider the total energy cost E(R, A<f>) = 
J G ff(A0) 2 /2 to create a small phase difference A<f) be- 
tween two faraway points A and B separated by distance 
R, where J e g is the disorder-averaged stiffness. For the 
classical XY chain with a random distribution of cou- 
plings P(J), it is known |2(J that there is a transition 
depending on the exponent P(J) oc J a as J — * 0. Fol- 
ey > 0, there is effectively a nonzero mean stiffness, and 
Jes(R) ~ R- 3/4 . For a < 0, there is a qualitatively 
weaker stiffness with exponent depending on a. 

Hence our model has a transition when a = £ , or 
when the typical length of a ID link on the percolation 
backbone is equal to the correlation length of the ID JJ. 
There are estimates of the number a for some standard 
lattices [l^. The stiffness transition is not at the same 
value of Ej / Ec as the KT transition where the ID cor- 
relation length becomes finite. The transition is caused 
by competition between the ID correlation length ( and 
the percolation physics that controls a. 

To improve the model, the 2D blobs can be modeled as 
2D classical percolation clusters rather than single mo- 
ments. (This is still an approximate model, since it ne- 
glects quantum effects in the blobs, of course.) With this 
modification, the asymptotic behavior of the sum is domi- 
nated by the blobs in the stiff phase: now E 7T (R) ~ R~ t l v 
in the stiff phase, where t « 1.31 is the classical perco- 
lation resistivity exponent and v = 4/3. The location of 
the transition to the floppy phase is also modified. The 
conclusion is that, at least in this approximate model, 
there is a stiff phase described by E 7r (R) ~ Rr l l v but 
with only ALRO, or even SRO with sufficiently large £, 
on the ID links. Quantum fluctuations on the links do re- 
duce the stiffness numerically in this stable "stiff" phase, 
and in the toy model can induce an exotic "floppy" phase. 

In the physical system where the blobs also are 
quantum-mechanical, the sharp transition between stiff 
and floppy phases will be smeared out since there is no 
obvious order parameter to protect it; and there will ap- 
pear a disordered phase for large Ec- This model may 
explain the numerical observation of || on the spin-half 
Heisenberg model that the stiffness has the percolation 
exponent but is numerically less than the classical value. 

The above description of the percolation cluster de- 
pended only on the assumption of either algebraic order 
or a correlation length in the one-dimensional links, so 
the same physical picture describes also the O(N) rotor 
models for N > 2 and the Heisenberg quantum antiferro- 
magnet. For the antiferromagnet, the ID links are alge- 
braically ordered for half-integer spin with short-ranged 
interactions, but develop a correlation length for either 
integer spin or strong frustrating interactions. 

The remainder of this paper discusses observable con- 
sequences of the picture suggested by the above model: 
that there is a stiff or "renormalized classical" phase 
where quantum effects modify the system quantitatively 
(by reducing the stiffness) but not qualitatively. Such a 
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FIG. 3: The three different stiff-phase forms for the correla- 
tion length A near the T = 0, p = p c geometric critical point. 
Above p c the correlation length is just set by the percolation 
correlation length (cluster size) £. Below p c , p a vanishes as 
\p-p c \ x , x » 1.31. 

model has been previously used to fit T = numerics at 
Pc 0- and similar renormalized classical physics applies 
to the undiluted 2D AF 0|. Here we extend the renor- 
malized classical picture to T > and p < p c (Fig. 0. 

At Pc, the correlation length A diverges as a power- 
law as T — > rather than as an exponential, and the 
power-law is given by the classical resistivity exponent: 
A ~ (J/kT)"/ t . The classical Heisenberg model (i.e., 
s a classical unit vector in Q) has a crossover from 
threshold behavior to quasi-2D behavior for p < p c . 

We can estimate this crossover using a one-step 
real-space RG. The length scale that characterizes the 
crossover from threshold to 2D physics is the percolation 
length £ ~ \p — p c \ v ■ Estimating the coupling strength 
at this length from the threshold stiffness, and using this 
coupling strength and length scale in the standard ex- 
pression for A in the renormalized classical 2D AF, gives 
for the correlation length 
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which is shown in Fig. Here, v = 4/3 and p s (p) ~ 
J\p — p c \ x , with x — t sa 1.31. This estimated value for x, 
which in our one-step RG controls the effective coupling 
in the 2D AF as p — ► p c , has previously been conjectured 
for the bulk modulus as p — > p c • 

There are several sets of experimental and numerical 
data to which the above forms can be compared. The 
T = 0, p = p c numerics of Sandvik do show classical 
percolation scaling, which is partly explained by our con- 
struction of a model stable to quantum fluctuations in the 
links. Correlation length measurements on the DHAF 
from neutron scattering Q seem to exhibit a crossover 
with dilution. The crossover was fit to flfl | 



e c/a 
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(12) 



a 8 2np s 1 + {4n Ps /T)-»T- 
At high-T, (|12|l like the classical model, follows a power- 



law rather than an exponential. An intermediate scaling 
regime controlled by a bilayer multicritical point |ll| has 
been proposed Q to explain deviation of observed scaling 
[vt ~ 0.7) from the classical value (vt ~ 1.02). 

In conclusion, superconducting or antiferromagnetic 
order can exist to p c in 2D but has low stiffness because 
of both cluster geometry and quantum fluctuations. This 
reduced p s results in a correlation length much shorter 
for T > than without dilution. For diluted La2Cu04, 
our analytic results explain and compare reasonably well 
with existing experiment and numerics. 
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